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Recent Monte Carlo simulations of a grafted semiflexible polymer in 1 + 1 dimensions have re- 
vealed a pronounced bimodal structure in the probability distribution of the transverse (bending) 
fluctuations of the free end, when the total contour length is of the order of the persistence length 
[G. Lattanzi et al, Phys. Rev E 69, 021801 (2004)]. In this paper, we show that the emergence of 
bimodality is related to a similar behavior observed when a random walker is driven in the trans- 
verse direction by a certain type of shear flow. We adapt an effective- medium argument, which 
was first introduced in the context of the sheared random- walk problem [E. Ben-Naim et al, Phys. 
Rev. A 45, 7207 (1992)], in order to obtain a simple analytic approximation of the probability 
distribution of the free-end fluctuations. We show that this approximation captures the bimodality 
and most of the qualitative features of the free-end fluctuations. We also predict that relaxing the 
local inextensibility constraint of the wormlike chain could lead to the disappearence of bimodality. 
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Semiflexible polymers have been the on focus of in- 
tense theoretical and experimental research activity in 
recent years. The main reasons are their relevance to 
biology, as many biologically important macromolecules 
fall in this category (including the building blocks of the 
cytoskeleton) Jj, i^j , and inherent challenges in their the- 
oretical study. A widely used minimal theoretical model 
of semiflexible polymers is the wormlike chain, a locally 
inextensible, fluctuating line with bending stiffness [3|. 
The main parameter in the description of a wormlike 
chain is the persistence length, Lp. It is defined as the 
correlation length of the tangent unit vector along the 
polymer contour, and it is proportional to the bending 
stifFnes, k: Lp = 2K/[kBT{d — 1)], where T is the tem- 
perature, and d is the dimensionality of the embedding 
space. The flexible limit of the wormlike chain, L ^ Lp 
(L being the total contour length), and to a lesser ex- 
tent the weakly-bending limit, L ^ Lp, have been well 
studied theoretically. In both of these extreme cases, the 
wormlike chain exhibits scaling behavior. In the inter- 
mediate regime where the total contour length is of the 
order of the persistence length, however, scaling ceases to 
exist and intriguing new phenomena emerge. One of the 
most striking features in the fluctuations of a polymer 
in this intermediate regime was revealed in recent Monte 
Carlo simulations of a grafted wormlike chain in a two- 
dimensional embedding space |^] . A grafted polymer has 
both the position and the orientation at one end fixed. 
If we look at the probability distribution of the trans- 
verse (bending) displacement of the free end, it starts 
as a delta function in the rigid- rod limit, which devel- 
ops into a Gaussian in the weakly-bending regime, and 
it becomes Gaussian again (of a different type) in the 



flexible regime. Surprisingly, in the intermediate region, 
the probability distribution is not a smooth interpola- 
tion between the two Gaussian limits, but it displays a 
pronounced bimodality. One can safely claim that this 
bimodality is the hallmark of semiflcxibility. 

We should point out that the emergence of bimodality 
(or multimodality) from an initially unimodal probabil- 
ity distribution upon the variation of a time-like variable 
is not peculiar to the semiflexible polymers. It has been 
shown that Levy flights in steeper than harmonic poten- 
tials exhibit a critical time beyond which an initially uni- 
modal distribution evolves into a bimodal terminal one 
5] . Depending on the steepness of the confining poten- 
tial, a trimodal transient may exist 0]. Another mani- 
festation of bimodality which is formally closer to that 
in semifiexible polymers occurs in diffusion-convection 
problems. If a neutral Brownian particle is carried by 
a power-law shear fiow in the transverse direction, then 
the probability distribution of transverse displacements 
gives rise to a terminal bimodality, depending on the ex- 
ponent of the flow profile |3] . A similar bimodality occurs 
when the convective fiow is odd and random (for a single 
realization of the randomness) 0- The most important 
difference between the above mentioned cases and the 
probability distribution of transverse fiuctuations in the 
grafted semiflexible polymer is that, in the former, bi- 
modality is terminal, whereas in the latter it is transient. 

If G[r{L),9{L)] is the probability distribution func- 
tion for a wormlikc chain having its free end at point 
r(L) = [x{L), y{L)] with a tangent vector making an an- 
gle 0{L) with respect to the clamping direction (x-axis), 
given that the other end is at the origin of the coordinate 
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system, it obeys the equation 0: 
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If we integrate out the spatial degrees of freedom {x and 
y), the resulting diffusion equation for the angle is simply 
the differential equation associated with a path integral 
of the Boltzmann weight of the bending energy . The 
"convective" terms express the local inextensibility of the 
wormlike chain pH l. 

If we integrate out the longitudinal fluctuations, x(L), 
we obtain: 
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Gy[yiL),9{L)]=0. (2) 



Integrating out the angle in Gy[y{L),6{L)] would give 
us the desired distribution which was observed in the 
simulation of Ref. 01 • Eq. 10) can be simplified and 
solved analytically (using Fourier transformations) in the 
weakly bending limit (i <C ip), where 9 I and sin 9 ~ 
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where (w) = 9/2. We have introduced the extra sym- 
bol V because we want to invoke the analogy with the 
diffusion-convection problem of Ref. 1^. According to 
this analogy, Eq. Q describes the motion of a particle 
in a two-dimensionsional space {9 and y) , with L becom- 
ing the time variable: it diffuses in the ^-direction and it 
gets carried by a sinusoidal flow in the y-direction. The 
flow is v{9,y) = {sm9)y. Note that, in contrast to the 
flows considered is Refs. and 0, this flow has a pe- 
riodic proflle. The weakly-bending limit of the wormlike 
chain corresponds to the case of linear shear flow. In that 
case, (v) = 9/2 is the average velocity of the y-coordinate 
in the interval [0, 9]. 

Before we implement the effective-medium approxima- 
tion, we should mention that attempts to perturbatively 
capture the onset of bimodality (from the stiff limit) have 
failed. Specifically, we had expanded sin 9 in Eq. |2Jl up 
to the quintic term, and we had treated the nonlinear 
terms to lowest order (one-loop) in perturbation theory. 
We had also used the analytic expressions for the variance 
and the kurtosis of the exact distribution in an Edge- 
worth expansion [l^ about the Gaussian, again without 
success. 

In the effective-medium approach, following Refs. 
and 1^, we hypothesize that the probability distribution 
of the transverse displacement and the slope of the free 
end maintains the form of Eq. Q for arbitrary bending 
stiffness (arbitrary L/Lp), but with the average "veloc- 
ity" (v) now being replaced by an "effective velocity," 




FIG. 1: The distribution funcion of bending displacements, 
P{y), from the analytic approximation and from the simula- 
tion, approaching bimodality from the the stiff limit, y has 
been rescaled to y/L and s = L/Lp. The lower panel shows 
the onset of multimodality with the three peaks. 
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FIG. 2: The distribution funcion of bending displacements, 
P{y), from the analytic approximation and from the simula- 
tion, as the polymer approaches the flexible region. The lower 
panel shows a detail from the region where bimodality is most 
pronounced. 



namely sin 9. In this approximation, the probability dis- 
tribution of transverse fluctuations reads: 
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The integral can be evaluated numerically, and the re- 
sult compared with the Monte Carlo simulation data is 
illustrated in Figs, ^andl^l In the weakly-bending limit. 
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the analytic approximation gives a narrow Gaussian with 
a width (y^) = (13/6)L'^/Lp (there is a discrepancy by 
a numerical prefactor of about 0.31 with the exact re- 
sult). As L/Lp increases, the distribution flattens and 
it develops three peaks at L/Lp « 0.77 in the analytic 
approximation and at L/Lp w 0.75 in the simulation. In- 
creasing the flexibility of the chain, the central peak gets 
suppresed, and a clear double-peak structure emerges, 
which is most pronounced around L/Lp k, 1.5 both in 
the analytic approximation and in the simulation. The 
results for the position of the peaks from the two ap- 
proaches are quite close. Eventually, the bimodality gets 
suppressed and a unimodal distribution reemerges in the 
flexible regime. This happens at L/Lp w 2.8 in the sim- 
ulation and at L/Lp w 3.8 in the analytic approxima- 
tion. The reentry into unimodality is accompanied by a 
transient trimodality in the simulation whereas this tri- 
modality does not appear in the analytic approximation. 
The major drawback of the effective-medium approxima- 
tion is that it fails to yield the right scaling for the width 
of the Gaussian distribution in the flexible regime. It 
reproduces [y"^) ^ L^^"^ of the weakly bending limit, 
instead of \/{y^ ~ L^/^ of the Gaussian chain. It also 
has Gaussian tails which spread into a region of extreme 
dispacements which is forbidden by the inextensibility 
constraint of the wormlike chain. 

We now consider the probability distribution of the 
longitudinal displacements of the free end. If we integrate 
out the transverse displacements, y{L), in Eq. we 
obtain: 
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Applying the effective-velocity approximation to this 
equation, and integrating out the orientational degree of 
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Plotting P{x) for various values of L/Lp, we obain a 
distribution which, apart from the Gaussian-fat tails, 
qualitatively agrees with the results from the simulation. 
Specifically, the value of x which corresponds to the peak 
of the distribution is pretty close to that obtained from 
the simulation over a wide range of stiffness (Table I). 
The agreement is better in the weakly-bending region. A 
remarkable feature of this distribution is that, although 
it always has only one peak, for 1.8 < L/Lp < 3.4, it 
exhibits two convex bumps. This structure was predicted 
by the analytic approximation and was subsequently con- 
firmed by the simulation. The presence of the two bumps 
for L/Lp = 2.8 is shown in Fig. |2| Note that this struc- 
ture is somehow reminiscent of the double-peaked struc- 
ture in the spatial density function of a (free) wormlike 
chain which was observed in Monte Carlo simulations by 
Dhar et al. in Ref. i in an overlapping region of L /Lp. 
Using Eq. © to calculate the average longitudinal posi- 
tion of the free end, (x), we obtain a particularly simple 
expression: 
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In Table II, we show how this approximate result com- 
pares with the value of (x) obtained from the simulation. 
The two results are very close in the weakly-bending 
regime, but they diverge very quickly as the wormlike 
chain enters the flexible regime. 





0.1 0.2 


0.3 0.4 0.5 0.75 1 1.5 2.5 2.8 


simulation a^max/L 
analytic a;niax/L 


0.992(1) 0.984(1) 
0.967 0.943 


0.977(1) 0.968(1) 0.960(1) 0.940(2) 0.922(2) 0.878(2) .798(4) 0.762(4) 
0.925 0.908 0.893 0.860 0.827 0.770 0.684 0.653 




2.5(1) 4.2(2) 


5.3(2) 6.2(2) 7.0(2) 8.5(3) 10.3(3) 12.3(3) 14.3(6) 14.3(6) 



TABLE I: Longitudinal position of the free end at the peak of the probability distribution, from the simulation, from the 
analytic approximation, and their relative difference, for various values of the flexibility, L/Lp. 

L/Lp I 0.1 0.2 0.3 0.4 0.5 0.75 1 1.5 2.5 2.8 

simulation {x/L) 0.953 0.910 0.869 0.829 0.792 0.711 0.639 0.525 0.373 0.340 

analytic {x/L) 0.951 0.905 0.861 0.819 0.779 0.687 0.607 0.472 0.287 0.247 

5(%) pa2 05 09 LO L6 3^4 5i) 10.1 23.1 27.4 

I 



TABLE II: Average longitudinalpo- 
sition of the free end, from the simu- 
lation, from the analytic approxima- 
tion, and their relative difference, for 
various values of the flexibility, L/Lp. 



Besides providing simple approximations for the prob- placements, the effective-medium approach can also be 
ability distributions of transverse and longitudinal dis- used in order to provide an approximation for the com- 
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FIG. 3: The distribution funcion of longitudinal displace- 
ments, P{x), from the analytic approximation and from the 
simulation, for L/Lp = 2.8, where it exhibits two convex 
bumps. X has been rescaled to x/L. 
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FIG. 4: Color density plot representing the two-dimensional 
probability distribution, P[x,y), of the free end, for L/Lp — 
2/3, obtained from the analytic approximation (Eq.JHJ). The 
longitudinal position, x, and the transverse position, y, have 
been rescaled to x/L and y/L, repsectively. 



plete distribution function, G{x,y,9): 
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Integrating out the angle, we obtain the probability dis- 
tribution function for the two-dimensional position of the 
free end, P{x,y). This distribution, for L/Lp = 2/3, 
is shown in the density plot of Fig. ^ Comparing it 
with Fig. 2a of Ref. one can see that the effective- 
velocity approximation qualitatively agrees with the sim- 
ulation, and the results for the position of the crest ob- 
tained from the the two approaches are even quantita- 
tively close. The quantitative agreement improves as the 
polymer gets stiffer. 

Out of mathematical curiosity (and not only!), one 
may want to explore a modified version of the diffusion- 



convection problem described by Eq. ^ . The new prob- 
lem contains an extra, dimensionless parameter e as the 
amplitude of the convective velocity: 
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While Eq. has only one scale, the parameter e in 
Eq. ^ causes a separation of scales between the diffu- 
sion process which acts on the time scale t = L and the 
convection process which acts on the time scale T = L/e. 
A rigorous multiscale-perturbation-theory Q analysis 
of this problem is deferred to future work. The imple- 
mentation of the simple effective-medium approach, how- 
ever, yields quite interesting results. We now consider the 
effective- velocity to be equal to (v) — es'mO, and we look 
at the dependence of P{y) on e. It turns out that, as e 
decreases, the width of the i/Lp-interval which exhibits 
bimodality shrinks, and, at a "critical" value ec ~ 0.66, 
P{y) becomes unimodal for any L/Lp. In the context 
of the wormlike chain, softening the relative strength of 
"convection" can be interpreted as relaxing the local in- 
extensibility constraint. We leave the disappearence of 
bimodality with the softening of the local extensibility of 
the polymer as a conjecture, which will be elucidated in 
a future work. 

Summarizing, in this paper, we used the analogy be- 
tween the wormlike chain and a diffusion-convection sys- 
tem, and we applied an effective-medium approach to 
analytically account for some unexpected features in the 
distribution of the semiflexible polymer conformations. 
It is remarkable that this very simple approximation has 
been so successful for three qualitatively different types 
of convective flow: random flow 0, power-law shear flow 
0, and - in our case - periodic flow. 

P.B. thanks R. Metzler for interesting discussions. 
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